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Abstract 

For quantum systems described by Schrodinger operators on the half-space W^~^ x M-^ 
the boundary force per unit area and unit energy is topologically quantised provided the 
Fermi energy lies in a gap of the bulk spectrum. Under this condition it is also equal to 
the integrated density of states at the Fermi energy. 

1 Introduction 

Consider a quantum system on the half space M*^"^ x M-''. One distinguishes between its 
behaviour for Xd << and Xd near speaking about the bulk behaviour in the first, and about 
the edge (or boundary) behaviour of the system in the second case. Bulk and edge behaviour 
are not independent but topologically quantised observables in the bulk related to topologically 
quantised observables at the edge. A famous example of this type is the quantum Hall effect in 
which the Hall conductivity can either be related to a current current correlation in the bulk 
or is simply the conductance of the edge current iPfSSl EMI ISKHOni IKHSn2l IHX;n2l IKSOabj In 
this letter we present another example relating the value of the integrated density of states on 
a gap (i.e. at energies lying in the gap) of the bulk spectrum to a force the boundary exhibits 
on the edge states. 

We discuss in the next section the underlying model and provide a proof of our claims in 
the case of one-dimensional periodic systems. The proof for the general case will be given 
elsewhere. 

2 The model 

In line with recent descriptions of aperiodic systems |Pa8n[ IB85| IB92[ IKSOSaj we describe 
the bulk behaviour of the system in the one particle approximation by a covariant family of 
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Schrodinger operators {H^}^^q, 
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all acting on L^(M"'), where A is a vector potential for an external constant magnetic field and 
n is a space of disorder configurations or of configurations which cannot be macroscopically 
distinguished. The second case is interesting for ordered systems, like quasi-crystals, in which 
Q can be the hull of a single point pattern describing set of average positions of the atoms 
and possibly decorated with information to distinguish between the kind of atoms, fl carries 
three structures. First, a metric topology in which it is compact, namely two configurations 
are deemed e-close if they agree up to an error of order e on a ^-neighbourhood around the 
origin G M'' (see |FHK02j for a precise formulation). Second, an action of the group of 
translations which we denote hj u x ■ u where x ■ u is the translate of the configuration 
a; e f2 which looks around x like uj around 0. Then a family of potentials {V^ji^en is called 
covariant if Vi^{x — y) = Vy.Lo{x) which implies that the corresponding family of Schrodinger 
operators {H^}^^^ is covariant w.r.t. magnetic translations. Third, fl comes with a translation 
invariant Borel-probabiblity measure P and all measurable quantities are averaged over Q with 
respect to this measure. This average has the meaning of a disorder average but should also 
be carried out if the configurations in Q describe ordered systems, because they are supposed 
to be macroscopically indistinguishable. The union IJwGf^ <7{H^) of all spectra a{H^) is called 
the bulk spectrum. 

2.1 Integrated density of states on gaps 

The integrated density of states at energy E oi a, single operator if^ from the family can be 
defined as the trace per unit volume r of the spectral projection Pe{H^) of onto the states 
up to energy E. If the probability measure P is ergodic then the trace per unit volume of 
Pe{Hu)) is P-a.s. constant over and 



Thus the integrated density of states should be thought of as an expectation value of the whole 
family of operators. The integrated density of states is constant on gaps in the bulk spectrum. 
It has been shown that equation has a C*-algebraic interpretation if E lies outside the bulk 
spectrum. In particular, the family {-PB(-ff(^)}a;en can be viewed as a projection in the natural 
C*-algebra associated with the configuration space VL (and the magnetic field) and that the l.h.s. 
of depends only on the homotopy class of Pe in this C*-algebra |B92j . This formulation 
makes clear that the value of the integrated density of states on a gap is topologically quantised, 
namely, first, it is stable under perturbations of the Schrodinger operator by covariant operators 
which do not lead to a closing the gap, and second, it lies in a countable subgroup of M, the gap- 
labelling group, which depends only on the topology of fi, its measure P, and the translation 
action, but not on the specific form of the potentials. In specific cases this group can be 
determined, e.g. for hulls of Delone sets of finite local complexity this group is the sub-group 
generated by the relative frequencies of the patterns appearing in the Delone set [BBCiUlBmiKPj . 
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2.2 Boundary force per unit area and unit energy in gaps 



To describe the behaviour near the edge we consider the same family of operators but restrict 
them to L^(R'^~^ x M-*) demanding Dirichlet boundary conditions at the boundary M*^"^ x {s}. 
We define the total boundary force to be minus the variation of the energy under a variation 
of the position s of the boundary. To describe this in one and the same Hilbert space, say the 
one corresponding to s = 0, we vary instead the position of the configuration in the opposite 
direction. Per unit area the boundary force exhibited on the states in some interval A is 
therefore 

/ dP{uj)f{PA{HMHu.)) , KHu) ■■= lim ^ , 

Jn ^d^o Xd oxd 

where f is the trace per unit area parallel to the boundary and we have for clarity denoted the 
restriction of to L^(M'^~^ x M-°) with Dirichlet boundary conditions with a hat. Note that 
we include the P-average directly into the definition of the boundary force. This is crucial for 
its topological quantisation and it is not true that f{PA{H^)5{Hi^)) is P-a.s. constant over VL 
for ergodic measures. The boundary force per unit area and energy at energy E is then 

|A| is the length of the interval which tends to in that limit. The limit exists if E lies in a 
gap of the bulk spectrum. Moreover, if E lies in a gap then Fh{E) is a topological quantity. In 
fact, it can be shown that 
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dP{u;)f{PA{H^)6{H^)) = 27rf((W:(A) - 1)5(W.(A))) (3) 
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where 



W^A) - 1 = (^e2-*(^-^») - l) Pa(^J , t=TX\^ i?o = min A 



Thus U^{A) is essentially the time evolution of the states of energy in A by the time which 
is the inverse of the width of A. Equation Q can be interpreted in a C*-algebraic context 
and shown to depend only on the homotopy class of the unitary Ui_j{A). Therefore, Fb{E) is 
topologically quantised in the same way as the integrated density of states. In fact, we can 
show using the tools of non-commutative topology of C*-algebras developed in |KSn3bj that, 
for energies E in gaps and ergodic P, 

mE)\ = im{E). (4) 

The proof of this result in the general case will be given in a separate publication. In the next 
section we give an elementary proof of the above equality for periodic one-dimensional systems. 



3 One-dimensional periodic systems 

We consider in this section the probably simplest case, in which we have a one- dimensional 
periodic configuration ujq, of period L, and fl is the set of its translates (no external magnetic 
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field). Then Q = {x-ljo\x ^ M} = M/LZ with (normalised) standard action of M by translation 
and Lebesgue measure. If we choose a differentiable periodic potential V then V^.^^^^x) = 
V{x + ^) defines a covariant family of potentials. For simplicity we write in place of V^.t^g- 
We need to combine results about the spectral theory of the family := —d^ + V^, ^ E [0, L), 
on three different spaces, see e.g. jDS88[ IB92j for background information. 

1) On L^(]R). Fixing G M and C, G [0,-Z^) one finds for each a, 6 G C a unique solution 
of {H^ - £^)\l/ = with initial condition ^(0) = a, \l/'(0) = b. is a function over M 
which is not normalisable. Since is a linear operator these unique solutions define a 

linear map — > C^: ^ ^''jo) ) '""^ ^'(^) ) ^'^^^^ associated matrix is called the 

monodromy matrix, we denote it by M{E,C,). The spectrum of on L^(]R) is the set 
G M| - 2 < Mu{E, + M22iE, < 2}. It is independent of ^ and hence equal to the 
bulk spectrum. It is convenient to call the closures of the connected components of the 
open set {E eR\-2< Mii{E,^) + M22(^,0 < 2} the bands. Then bands may touch. 

2) On L^([0,L]) with Dirichlet boundary conditions at the boundary points. On that 
space the spectrum of is a discrete countably infinite set /^2(C)) ' " ' }; < 
/^j+i(^)). We call these spectral values Dirichlet eigenvalues of H^. They are determined 
by the equation M^i^E,^) = 0. Between and /ij+i(^) lies the j + 1th band of the 
bulk spectrum (counted from the lowest band). We call a Dirichlet eigenvalue /in(0 ^ l^ft 
or right eigenvalue if its corresponding eigenfunction ipn^s^ satisfies |'?/'^_^(L)| < iV'n.glO)! 
\ip'^^^{L)\ > |V'n,g(0)l; respectively. The terminology comes from the exponential increase 
if the functions are considered over many periods which physically means that eigenfunc- 
tions of right eigenvalues are localised at the right edge. If two bands touch they touch 
at a Dirichlet eigenvalue which is neither a left nor a right eigenvalue. 

Sturm-Liouville theory gives us the important information that a real eigenfunction ipn,£, 
of to yUn(0 has exactly n zeroes on the half open interval [0, L) (so the zero at L is 
not counted). 

3) On L^(M-°) with Dirichlet boundary condition at the boundary 0. The spectrum of 
on that space is the union of the bulk spectrum with the right Dirichlet eigenvalues. In 

fact, E belongs to that spectrum iff the eigenvalues of M{E, ^) have modulus 1 or ^ 

is an eigenvector of M{E,C,) to an eigenvalue p of modulus strictly larger than 1. In the 
first case E belongs to a band and in the second the corresponding eigenfunction satisfies 

K,cWI=pK,c(o)l>l<,e(o)|. 

We give an alternative description of left or right Dirichlet eigenvalues, namely finiO is a 
left or right eigenvalue iff /i^(^) > or /i^(^) < 0, respectively. In fact, we calculate /i^(0 
for right eigenvalues 

l^n{0=9^ dxijjn,^{x)Hi:ljjn,^{x) = / dx\lijn,^{x)\'^Vl{x) (5) 



where i/^n,^ is the normafoec? eigenfunction, J° dx\'4^n,^ix)\'^ = 1. Clearly, this normali- 
sation is only possible since /in(0 is a right eigenvalue. Using integration by parts and 
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i'n,iiO) = 'ipn,tii-oo) = ip'^^^i-oo) = WQ find, 

fO 



dx\ipn,ii{x)\ V^ix) = - dx4j'^^^{x)ilJn,i{x)V^{x) - c.c. 

J — oo 



^ ^ ^ 

dxij'„^{x){EiJn,^{x) + 7^i^n,^{x)) - C.C. 

oo 



For left eigenvalues one proceeds similarly, but uses the space L^(M-°) instead. Since is 
then the left boundary of the integral one obtains a relative minus sign in the calculation. 
The remaining case is that /in(0 is neither a left nor a right eigenvalue. Then it must be 
at a band edge and therefore an extremum of fin- 

This has the following consequence which is crucial below: Since a real eigenfunction to a 
right Dirichlet eigenvalue finiO has n zeroes on [0, L) the equations finiO = A*; /^n(0 < 
have exactly n solutions on ^ G [0, L). 

Using an approximation of the trace per unit volume by the trace per unit volume in the 
representations L'^{[—NL, NL]) with Dirichlet boundary conditions one obtains from Sturm- 
Liouville theory in the limit N —>■ oo 

IDS(/.„(0) = (6) 

independent of ^. 

We determine the boundary force If A is an interval in the rath gap of the bulk 

spectrum then has exactly one non-degenerate eigenvalue (a right-Dirichlet eigenvalue) 
provided /i(^) G A and < 0, otherwise it has none. Hence the integral kernel {x\P/^{H^)\y) 

of Pa{H^) is 4'n,i{x)4'n,i{y)XAif^n{0)^i~f^niO) whcrc Xa characteristic function on the 

interval A and G the Heavyside function. In one dimension f is the operator trace and therefore 

f(PA(4)n')= r rf^l^n,5(a;)l'XA(/^n(e))e(-/i;(0)n'(a:). 
J — oo 

By ^ and the fact that finiO = /^n(0 < has exactly n solutions on ^ G [0, L) 

r dxl^PnMW^if^niOM-fl'niml = r^^/^n(OXA(/^n(O)0(-/in(O) 



= — n dfj, = — n|A|. 
J A 

As it should be, this expression is negative and so the force points into the sample. We conclude 
that 

and so the strength of boundary force per unit area and unit energy is equal to the integrated 
density of states for an energy E which lies in a gap of the bulk spectrum. We end this section 
with some remarks. 
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1. The integer n appearing in the above expression may also be interpreted as a winding 
number of the Dirichlet eigenvalue on the complex spectral curve of H^. This is similar 
but not equal to the phenomenon observed in |H93j . 

2. The above also yields, for arbitrary interval A in the nth gap, 

Here Wn is the width of the rath gap and, for the +-sign ( — sign), ipn,£, is a normalised 
eigenfunction to the left (right) Dirichlet eigenvalue. This equation seems interesting in 
its own right for one-dimensional periodic operators. 

3. Since, by the boundary conditions. 



x=0 



we see that the boundary force per unit area and unit energy is determined by the P- 
average of the first non- vanishing coefficient in the Taylor expansion of the density of the 
particles at the edge. 

4. We note again that the P-average is crucial for the topological quantisation. In fact, 
^n,?(0) tends to if ^ tends to the extrema of the function /Xn(0- 
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